Surface influence on stability and structure of III-V nanorods: First-principles studies 
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We report ab initio investigations of hexagon- shaped, [111]/[0001] oriented III-V semiconductor nanowires 
with varying crystal structure, surface passivation, surface orientation, and diameter. Their stability is dominated 
by the free surface energies of the corresponding facets which differ only weakly from those of free surfaces. 
We observe a phase transition between local zinc-blende and wurtzite geometry versus preparation conditions 
of the surfaces, which is accompanied by a change in the facet orientation. The influence of the actual III-V 
compound remains small. The atomic relaxation of nanowires gives rise to smaller bond lengths in comparison 
to the bulk znic-blende structures resulting also in somewhat reduced bilayer thicknesses parallel to the growth 
direction. 

PACS numbers: 61.50.Ah, 68.35.Md, 68.47 .Fg, 68.70,+w 



I. INTRODUCTION 

Highly anisotropic needlelike crystals (whiskers) have long 
been subject of physics and materials science. This inter- 
est, especially in the ultimately thin varieties, has been re- 
cently stimulated by the potential need as building blocks 
for nanoscale electronic and photonic devices 1,2 . The syn- 
thesis of nanowires based on a wide range of material 
system s 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 have been reported. Even heterostruc- 
tures have been prepared as vertically segmented wires along 
the growth axis or as core-shell structures perpendicular to 
the wire orientation. 1112 Diameters varying from about one 
to hundred nanometers and lengths extending to several mi- 
crometers are observed. Because of their considerable tech- 
nological potential for optoelectronics or high-speed electron- 
ics, nanowires or nanorods consisting of III-V semiconduc- 
tors, AlAs, GaAs, InAs, GaP, and InP, are of special interest. 
In the most cases the growth direction of III-V semiconductor 
nanorods is parallel to a [1 1 1] axis of the bulk zinc -blende (zb) 
substrate, which is the ground-state structure of the common 
III-V compounds. 

However, in contrast to bulk crystals the crystal struc- 
ture of the nanowires may differ noticable, depending on 
growth conditions and growth method. In particular, changes 
of the crystal symmetry from cubic to hexagonal stacking 
of the cation-anion bilayers have been observed in many 
case a 3 ' 4 ' 7 ' 8 ' 9 ' 10 ' 12 . There are reports of a rotational twin struc- 
ture around the [111] growth axis or even of the wurtzite (w) 
structure with [0001] stacking directio n 3 ' 4 ' 12 (schematically 
indicated in Fig. 1 ). This has been discussed in detail for GaAs 
and the unique vapor-liquid-solid growth mechanism— The 
growth temperature has an strong influence in the case of the 
metal-organic vapor-phase epitaxy (MOVPE). Free-standing 
GaAs nanorods possess zb structure for the temperature range 
460-500°C but changes to the w polytype at 420°C or above 
500°C. 2 In the case of the chemical beam epitaxy (CBE) the 
w structure has been observed at 540°C or below 400°C. 13 
Laser-assisted catalytic growth at higher growth temperatures 
usually results in the zb phased 

The consequences of the stacking variation for the geom- 
etry parameters such as bond lengths are not clarified. No 
significant changes have been reported for the bilayer thick- 



passivating atoms 




FIG. 1: (Color online) Schematic picture of a hexagon-shaped 
nanorod with diameter D and a certain length L, which is a multi- 
ple of irreducible slabs. 



ness in [111] growth direction. 5610 Only a very small reduc- 
tion of the distance between (111) lattice planes compared to 
bulk values was observed in transmission electron microscopy 
(TEM) studies for GaAs 5 ' 6 and InP. 10 Other studies report a 
slight increase of the c/a ratio discussing InAs in terms of a 
hexagonal geometry,^ 5 - 

Another not well clarified point concerns the favorable ori- 
entation of the surface facets. There seems to be agree- 
ment about the hexagon shape of the rod cross-sections with 
equivalent facet orientation ! 3 ' 7 ' 8 ' 10 ' 14 ' 16 For GaAs rods few pa- 
pers report the orientation of the side facets, either {112} 
(zb) I {1100} (w) or {110} (zb) I {1120} (w)M Scanning 
electron microscopy investigations found hexagonal cross- 
sections with {112} side facets for GaP nanorods with zb 
stacking, 17 while {110} side facets have been reported for 
InAs and InP nanorods . 10 ' 14 ' 16 However, two-dimensional de- 
fects such as twin segments or altering orientation may lead to 
wires that are microfaceted. 18 Growth of wires with triangu- 
lar cross-sections may also occur as found for GaAs rods with 
{112} side facets grown in [1 1 1]A direction. 17 Relatively lit- 
tle is known about the surface bonding and the contamination 
of the side facets. Because surface oxides such as In203 on 
InAs desorb above 500°C, one thinks that annealing at lower 
temperatures results in incomplete desorption of the surface 
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contamination. 4 Of course, growth in an ultrahigh vacuum 
system should result in almost contamination-free surfaces. 
However, also amorphous overlayers have been observed and 
attributed to oxide layers that form after GaAs is exposed to 

The growth of quasi-one-dimensional III-V semiconductor 
rods with varying crystal structure and bilayer distances raises 
the question of the limit of a bulklike description of bonding in 
these nanowires. Especially for small diameters the predomi- 
nance of surface atoms over inner, bulklike atoms will eventu- 
ally lead to bonds and, hence, geometries different from those 
of the bulk systems. Thereby, the details of such changes may 
be influenced by the reconstruction/relaxation of the clean sur- 
faces or the stoichiometry/thickness of the surface contamina- 
tion. Despite varying numbers of dangling bonds and a possi- 
ble passivation of these dangling bonds, trials to understand 
the surface influence are almost missing. The interplay of 
bonding arrangement and stoichiometry of rod surfaces and 
zb — w phase transition has not been studied so far. 

On reason for this lack is the limitation of first princi- 
ples calculations to rather small (thin) systems. This ham- 
pers the comparability with experimental results. Parameter- 
free calculations were performed mainly on carbon or silicon 
nanowires/tubes. The investigated structures had diameters up 
to about one nanomete r 19 ' 20 , which is close to the thinest ex- 
perimentally observed structures. However, to make reliable 
predictions for the zb — w phase transition in III-V nanowires 
it is necessary to deal with larger systems. 

In the present article, we try to overcome the limitations of 
ab initio calculations by combining two approaches, micro- 
scopic total-energy studies of ultrathin nanorods with vary- 
ing surface facet orientations, bonding and overlayers with 
macroscopic (thermodynamic) studies based on surface ener- 
gies which are computed for individual surfaces. We demon- 
strate rapid convergence of the rod formation energy per sur- 
face area of III-V semiconductor nanowires with rising diam- 
eter. This fact enables us to extrapolate our results to diame- 
ters observed experimentally. In particular, we examine free- 
standing, hexagon-shaped GaAs, InP, and InAs nanowires 
with varying surface facet orientation and passivation, several 
diameters, and two crystal structures close to the thermal equi- 
librium. The results are extrapolated to the case of extremely 
thick rods with almost bulklike surfaces. The favored geome- 
tries are studied versus the preparation conditions of the wires. 



II. MODELLING OF NANORODS 

A. Total energy calculation 

The total-energy (TE) calculations are performed in the 
framework of the density functional theory (DFT) and the lo- 
cal density approximation (LDA) as implemented in the Vi- 
enna ab initio simulation package (VASP) 21,22 . The outer- 
most s-, p-, and (in the case of In) af-electrons are treated as 
valence electrons whose interaction with the remaining ions is 
modeled by pseudopotentials generated within the projector- 
augmented wave (PAW) method 23 . The energy cutoff of the 



plane-wave basis is tested to be sufficient with 15 Ry. To 
model nanowires we use the common supercell approach 24 . 
Neighboring rods are separated by a vacuum region of about 
2 nm to minimize artificial interactions across the periodic 
boundaries. The Brillouin-zone summations are restricted to 
MxNx2 meshes of special points according to Monkhorst 
and Pack 25 , with M=N=5 for diameters D (as defined in 
Fig. 1 as the vertical distance of two parallel surface facets) 
smaller than 1 nm and M=N=3 for diameters D larger than 
1 nm. In the case of bulk zinc -blende structures the minimiza- 
tion of the total energy leads to theoretical cubic lattice con- 
stants a = 5.61 A (GaAs), 5.83 A (InP), and 6.03 A (InAs) 
as well as to negative cohesive energies per cation-anion pair 
(without spin-polarisation corrections) /x]}" lk l/ = —9.68 eV 
(GaAs), -9.49 eV (InP), and -8.88 eV (InAs). 

Free clean and passivated surfaces of zb or w crystals are 
studied within the repeated-slab method using material slabs 
of about 14 to 17 A separated by a vacuum region of the same 
thickness 26 . The material slabs usually represent multiples 
of irreducible slabs. For the {1 !0}zb surface we use four ir- 
reducible slabs each of them consisting of two atomic lay- 
ers. The {112}zb surface is described using three irreducible 
slabs with six atomic layers. All considered surfaces are stoi- 
chiometric and electrostatic neutral. Nevertheless, in the case 
of {112} facets the different numbers of dangling bonds per 
anion or cation give rise to two unequal surfaces. In these 
cases we can only calculate an averaged value for the two 
surface energies. Their application to the nanorod energetics 
does, however, not influence the results since such inequiva- 
lent surfaces occur only pairwise in form of opposite facets on 
a nanorod. To model the wurtzite surfaces in the same manner 
we apply 1 1 atomic layers for the { 1 100} w surfaces and seven 
atomic layers for the {1 120}w surfaces in the slab. All calcu- 
lations are performed using a lxl translational symmetry. A 
possible surface reconstruction is not taken into account. 



B. Nanorod structure 

The one-dimensional translational symmetry in growth di- 
rection [111 ]/[0001 ] is given by irreducible slabs of six atomic 
layers (three cation-anion bilayers) of zb in [111] direction 
or of four atomic layers (two cation-anion bilayers) of w in 
[0001] direction 26 . Therefore, sometimes one speaks about 
the 3C(zb) and 2H(w) polytype of a compound. In our sim- 
ulations we use slabs of the length L = 2\/3flo of two (three) 
irreducible zb (w) slabs. The diameter of the rod is charac- 
terized by the distance D between opposite surface facets (see 
Fig.l). Nanorods with hexagonal cross-sections and two dif- 
ferent stackings (zb or w) of bilayers in growth direction are 
studied as observed in most of the experimental studies of 
III-V semiconductor a 8 i 1 °i 16 , Therefore we restrict ourselves 
to those rods, where the point-group symmetry C^ v is con- 
served for arrangements of both types of stackings as in the 
zb- and w-structures. We do not consider nanowires with 
high-index side facets, because of their roughness and pos- 
sible higher dangling-bond densities. Surface normal orien- 
tations are studied only parallel to high symmetry directions. 
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FIG. 2: (Color online) Schematic stick and ball representation of 
adjacent facets on type-I (left-hand side) and type-II (right-hand side) 
nanowires. 



Hence, we end up with two different types I and II of side 
facets. Type-I wires contain {112} and { 1 100} side facets in 
the zb and w case, respectively (Fig. 2, left). GaAs rods with 
such surfaces can, e.g., be grown by MOVPE 3 . In the type-II 
case, { 1 10} facets for zb and { 1 120} facets for w occur as also 
observed experimentally (Fig. 2, right) . 7 ' 10 ' 14 

To construct the nanowires on an atomic scale we start from 
a cation-anion bond parallel to the growth direction in the 
center of the nanorod. The nearest-neighbor atoms occur at 
the corners of a tetrahedron along cubic body diagonals. In 
the wurtzite case each second tetrahedron in stacking direc- 
tion is twisted by 180° resulting in a shortening of the pe- 
riod of translational symmetry in [0001] growth direction to 
two cation-anion bilayers compared to the three bilayers in 
the zinc -blende case. 

The side facets (top view) of the two wire types considered 
are also shown in the lower panel of Fig. 3. They are stoi- 
chiometric and electrostatically neutral. Even the numbers of 
cation and anion dangling bonds (DBs) are equal in each irre- 
ducible slab. The absolute DB densities of the corresponding 
infinite surface facets (i.e., neglecting the edge influence) vary 
between for type-II zb- and w-wires, 4v / 2/ [«o v^] for 

type-I zfo-wires, and 3\/6/[2flo] for type-I w-wires assuming 
the same bonding tetrahedra for w stacking with ideal hexag- 
onal lattice constants a — ao/V2 and c — 2oq /V3. In the latter 
case the DB densities vary with the termination of the surface. 
Since the considered rods always exhibit both types of ter- 
minations at opposite facets we have just given the averaged 
value of both DB densities. To simulate a possible passivation 
of the DBs at the wire surfaces due to adsorption of environ- 
mental atoms or metal atom excess as schematically shown 
in Fig. 1, we use fractionally charged pseudohydrogen atoms 



with a valence charge Z = 0.75 (H*) or 1.25 (H**) 2 !. To ob- 
tain the equilibrium positions of the passivating pseudoatoms 
we apply a conjugate gradient algorithm based on first princi- 
ples forces. 

Figures 2 and 3 indicate remarkably different geometries of 
the two types of side facets. In the type-II case the side facets 
have only one DB per surface atom and exhibit rather sim- 
ilar geometries. For {110}z£> facets the cation-anion bonds 
form a staircase-like structure, while the twist of the bonding 
tetrahedra in the {2110}w case gives rise to zigzag chains in 
[0001] direction. In contrast, type-I side facets differ remark- 
able in their bonding structures. In the average, each surface 
atom possesses 4/3 (zb) or 3/2 (w) DBs. However, already for 
zb stacking their distribution depends on the orientation of a 
{112} facet. Neighboring facets, for instance (112) and (211), 
are different with respect to their DB distribution as indicated 
in Fig. 2. On the surface only cation-anion bonds parallel to 
[0001]/[1 11] occur. The atoms in the uppermost atomic layer 
possess an averaged value of 1.5 DBs for both crystal struc- 
tures. However, the atoms in the second atomic layer of the 
(1 12)zb surface exhibit one additional dangling bond. In prin- 
ciple, the situation in nanorods is somewhat more complicated 
because of the existence of edges which modify the picture of 
free surfaces for finite wire sizes. However, for thick rods, i.e., 
in the limit D — > °°, the facets represent free surfaces. 

C. Formation energy of a nanorod 

In order to study the stability of the different nanowires we 
investigate the formation energy of a nanorod E*°i with a cer- 
tain surface area consisting of symmetry-equivalent facets de- 
fined by diameter D and length L (defined in Fig. 1)^ Taking 
into account that in a rod not all group-Ill and group-V atoms 
are fourfold coordinated but the cation and anion DBs occur 
pairwise the formation energy obeys the expression 

- N mol v™ l - H «-N moi ALiH*-H», (1) 

where E^^Nm-v ,N mo i) is the total energy of a rod contain- 
ing Nui-v cation-anion pairs with DBs or passivated by 7V mo i 
pairs of pseudohydrogens H* and H**. This energy is related 
to the corresponding bulk energy Nui-v Mm-y an< ^ to trie en ~ 
ergy of the 7V mo i H* — H** molecules used to passivate the total 
number of 27V mo i DBs. We assume that the rods are in thermo- 
dynamic equilibrium with a zb III-V crystal that represents a 
reservoir for group-Ill and group-V atoms. The passivation is 
governed by a reservoir of H* — H** molecules with the chem- 
ical potential /!//*-#**■ The variation of this chemical poten- 
tial allows to model the environmental influence during the 
growth or annealing process as have been demonstrated for 
individual surfaces. 28 Here (calculated as the neg- 

ative molecule binding energy) corresponds to the situation 
where the rod surface is exposed to molecular pseudohydro- 
gen at vanishing temperature. However, due to the need to 
overcome the dissociation barrier, the stable passivated sur- 
face structures will not form immediately. The surface cov- 
erage depends on the environmental conditions or the partial 
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FIG. 3: (Color online) Schematic stick and ball representation of type-I (left-hand side) and type-II (right-hand side) nanowire models. In 
the case of the top views (upper panels) only the atoms in the irreducible slab in different atomic layers are shown. The top views of the 
corresponding surface facets (lower panel) represent one possible arrangement versus the length of two (three) irreducible zb (w) slabs. The 
circle size indicates atoms in different atomic layers 1-6, while the circle color represents cations (A) or anions (B). Pseudohydrogen atoms 
are represented by very small circles. 



pressures and temperatures in the regions. This dependence 
A)1h*-h™ of Hh*-h** can be approximated similarly to that 
of a two-atomic ideal gas. A slightly negative potential (e.g. 
Ajj,H*-H™ ~ — 1 eV in the case of true hydrogen) may model 
typical MOVPE growth conditions while fully passivated sur- 
faces may form for A/X//*_#« > where atomic pseudohydro- 
gen is available^ In a real preparation process, e.g. for oxi- 
dation of rod surfaces in air, the adsorption processes are more 
complex. In addition modifications of the surface stoichiom- 
etry may appear. Nevertheless, the study of the rod formation 
energy (1), which is based on Kramer's grand canonical po- 
tential, should give important information about the equilib- 
rium state of rods in a certain gaseous environment. A gener- 
alization of (1) towards individual group-Ill or group- V atoms 
is possible by introducing the corresponding chemical poten- 
tials, in order to model different ratios of these species during 
the growth process i 26,29 ' 30 

The formation energy (1) depends on the geometry (di- 
ameter D, length L) of the considered rod and the prepara- 
tion conditions, i.e., Aji H *_ H **. We refer this energy to the 
surface area A roc j = 2\/3 DL of the considered rod and set 
Ajj,H*-H** = 0. The resulting surface free energy per unit 

3132 

are a ' 

?N d = [Eg?(Nm-v,N moi )-N m -vliiii-v 

-WmolM™V]Mrod (2) 

is influenced by the energetics of the inner rod atoms as will 
be discussed later, the energetics on a given surface facet, and 
the energetics of the edges between two facets. It still depends 
on the thickness of the rod, which determines the relative con- 



tribution of the three energies. Instead by D, in (2) we charac- 
terize the thickness by the dimensionless parameter Nj : 

N * = \ \— (type I rods), (3) 
V 3 a 

N d = V2— (type II rods), (4) 

in terms of the ideal lateral lattice constant an/ V2. The factor 
2/y/3 is related to the 30° rotation between both types of rods. 

III. RESULTS AND DISCUSSION 

A. Energetics and phase diagram 

For Nd — > °° the energy (2) converges against the surface 
energy y = % of the clean or passivated surface with consid- 
ered orientation as demonstrated in Table I. Reconstruction or 
relaxation of the surfaces are not taken into account.— The 
convergence is faster for passivated rods than for rods with 
clean facets. For the w-crystal structure y^. (2) is not exactly 
equal to the free surface energy. It contains contributions re- 
sulting from the change of the crystal structure. However, 
because of the small difference in the cohesive energies, e.g. 

A l /// lk v( vv ) ~ J u /// lk v ( zb ) ~ 25 meV for InAs > th ose contribu- 
tions are very small and can therefore be neglected (at least 
for the considered numbers Njji-v)- Our surface results are in 
good agreement with other first principles calculations. For 
instance the free surface energy of the relaxed cleavage (110) 
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surface of InAs with the value of 40.5 meV/A 2 agrees well 
with the result of Ref.[H 

Table I exhibits three clear trends. First, the surface ener- 
gies of the rods are only slightly larger than the one of the cor- 
responding free surface. With rising rod diameter, N d — ► °°, 
a rapid convergence to the free surface values is observed. 
Already for InAs nanorods with diameters of about 1.7 nm 
(A^/=4, type II) the results are only slightly above those for 
very thick rods. Hence, the diameter of the rod is of little in- 
fluence on the stability of the rod, taking as measure the min- 
imum surface energy with respect to the surface area A ro d for 
given volume 33 . Second, for the assumed preparation condi- 
tions, Ahh*-h** — 0, the energies of the passivated surfaces 
are considerably smaller, indicating the stability of the passi- 
vated rods versus such with clean surfaces. Third, the type-I 
rods with a higher dangling-bond density of the surface pos- 
sess a larger surface energy independent of the rod diame- 
ter and the passivation. All values in Table I follow the se- 
quence of DB densities: type II(zb)«type II(vv) < type l(zb) 
< type I(w). A clear favorization of a certain stacking with 
a certain diameter, as recently reported by means of total en- 
ergy calculations using empirical pseudopotentials, 35 cannot 
be concluded. Which stacking order, zb or w, in the rod is 
more stable depends on the passivation and the wire type. 



a strong influence of the preparation conditions on the rod 
stability. If mainly passivating molecules of a not too high 
density are present, the rods prefer clean surfaces of type II 
with low dangling-bond density, i.e., { 1 10} / {2110} facets. 
In the presence of dissociated molecules (free radicals H* and 
H**), i.e., S.\Iu*-h** > 1 eV , the rod should be passivated 
and possesses type-I facets, i.e., those with originally higher 
dangling-bond densities. In this case, a clear favorization of 
the wurtzite stacking versus the zinc -blende stacking can be 
stated. However, in all other cases the energies of the two dif- 
ferent bilayer stackings are very close, so that probably other 
(e.g. kinetic) aspects of the rod growth are more important for 
a given stacking than energies derived for pure equilibrium 
situations. Under intermediate conditions, i.e., Ahh*-h** be- 
tween -1 eV and zero, passivated type-II facets are preferred. 
Although type-II facets have the same dangling-bond density 
for zb and w stacking, we find a slight preference of the zb 
structure. This is a consequence of the higher negative bulk 
cohesive energy $f} k v (w) of the w structure and the influ- 
ence of edges on the DBs. The experimental results concern- 
ing InAs are usually not unique. 4 6 However, in Refs.RlandfT^ 
both types of stacking, zb and w, have been observed under 
MOVPE conditions, which is in agreement with our predic- 
tions. 



TABLE I: Surface energies y^ d (in meV/A 2 ) for InAs rods of type I 
or II for varying diameter ~ Nj . The values for Nd = °° are calculated 
for free surfaces using the repeated-slab method. The first (second) 
value has been obtained for clean (passivated) nanorods. 
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The stability of InAs rods with different stacking, sur- 
face orientation, and surface passivation as a function of 
the preparation conditions is demonstrated by the phase di- 
agram in Fig. 4. According to a Wulff-like approac h 33,34 
the most stable rod should have the lowest surface energy 
YN d = ^surfMrod- This allows to identify the most favor- 
able facet orientation, surface coverage, and bilayer stack- 
ing for given preparation conditions and fixed diameter. We 
have plotted the surface energy YN d versus the variation of 
the chemical potential of the passivating species. Negative 
values of AjiH*-H** stand for passivating species-poor con- 
ditions, while positive values indicate the availability of free 
passivating atoms. Because of the weak diameter dependence 
of the energies jN d we use the free surface values /«, scaled 
with the corresponding DB densities, which are diameter de- 
pendent. In this way we are able to extrapolated our results to 
rods with a realistic diameter of about 5 nm. One observes 



TABLE II: Surface energies f^ d (in meV/A 2 ) for passivated GaAs, 
InP, and InAs rods of type I/II for varying diameter ~ N^. The values 
for Nd = °° are calculated for free surfaces using the repeated-slab 
method. 



N d 


c 

zb 


3aAs 

w 


zb 


InP 

w 


1 

zb 


nAs 

w 


2 

4/3 

OO 


25/21 
23/19 
21/16 


30/22 
28/22 
29/19 


26/22 
24/20 
24/17 


30/21 
28/22 
30/19 


31/28 
29/26 
27/23 


37/27 
33/28 
33/25 



In order to study chemical trends within the III-V semicon- 
ductors GaAs, InP, and InAs with zb bulk ground states, we 
restrict the computations to passivated rods according to the 
wide stability range in the InAs phase diagram in Fig. 4. The 
energy variation with the compound may be larger for clean 
surfaces. Nevertheless, they should only occur under ultra- 
high vacuum conditions. Results for % d are listed in Table 
II for different wire thicknesses. The least(most) stable pas- 
sivated nanorods are obtained for InAs(GaAs). However, the 
most important consequence of Table II is, that the surface 
energies of all considered semiconductors show the same be- 
havior and only a small reduction from InAs via InP to GaAs. 
Therfore the phase diagrams of the considered III-V semi- 
conductor nanorods will have the same qualitative shape as 
that of InAs in Fig. 4 Hence, the occurrence of both stack- 
ing types, wurtzite and zinc blende, in III-V semiconductor 
nanorods may be explained by the weak dependence of the 
energy sequence on the bilayer stacking and the rod diameter 
(as discussed in detail for InAs). 
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FIG. 4: (Color online) Surface energy per surface area y^ d of a InAs 
rod with a thickness of about 5 nm versus the variation of the chem- 
ical potential hfln'-H" of the passivating molecules. Solid lines: 
zb stacking, dashed lines: w stacking. Rod type and surface cover- 
age are indicated. The central part surrounded by thin dotted lines 
describes the region of the most reliable chemical potentials. 



B. Geometry and bond lengths 

All results presented until now have been computed without 
lattice relaxation. The atomic positions were obtained from 
the tetrahedral structure and the (theoretical) bulk cubic lat- 
tice constant ao- Only the stacking has been changed. In or- 
der to study the influence of atomic rearrangements due to the 
surface presence, we allow atomic relaxation in the nanorods 
under the constraint of conservation of the tetrahedral coor- 
dination. Since both stackings, zb (3C) and w (2H), can be 
described in hexagonal unit cells with lattice constants a and 
c, we allow for changes of these lattice constants with respect 
to the corresponding ideal bulk values a = ao/V2 (zb, w), 
c = \/3ao (zb) or c — -^ao (w)- For thin rods (A^=2) with 
clean side facets we find a nearly vanishing increase of the a 
lattice constant by 0.06 % for both stackings. Contrariwise the 
c lattice constant is reduced by 1.7 % (zb) or 2.8 % (w). As 
thermodynamically expected this leads to a small shrinkage of 
the III-V bond lengths in the rods. Surface tension is accom- 
panied by a tendency for volume reduction. A similar ten- 
dency has been observed for semiconductor nanocrystals. 36 
The absolute surface energies YN d are somewhat changed with 
respect to the unrelaxed rod, but their relative differences re- 
main almost the same. With the c lattice constant the cation- 
anion bilayer thickness parallel to the growth direction is de- 
creased. 

For passivated rods the already small changes nearly van- 
ish, i.e., the a lattice constant remains constant, while the 
c lattice constant is reduced by less than 0.8 % for both 
stacking types. This observation is in good agreement with 
the almost bulk-like bilayer thicknesses measured by high- 
resolution TEM 5 ' 6 ' 10 . For zb GaAs nanorods only a tiny re- 



duction (less than the error bar) of the lattice plane spacing 
in growth direction was found in Refs. 150, while for w-InP 
nanorods a reduction of the c lattice constant by 2.8% was 
observed. However, we cannot confirm the increase of the 
c/a ratio with respect to the ideal value, which was found in 
Refs.fliandfil 



IV. SUMMARY AND CONCLUSIONS 

Using ab initio total-energy calculations for thin nanorods 
and individual surfaces of the same materials we have stud- 
ied the energetics of III-V semiconductor nanorods in ther- 
modynamic equilibrium. For rod diameters smaller than two 
nm formation energies are directly calculated for a given seg- 
ment length. For thick rods the rod surface energy is obtained 
from values of the individual facets using a Wulff-like con- 
struction. The preparation conditions, especially the presence 
of passivating species, are simulated by corresponding chem- 
ical potentials (changing from the free energy to the grand 
canonical thermodynamic potential). We studied rods grown 
in [1 1 1]/[0001] direction with varying diameter, bilayer stack- 
ing, surface coverage, and facet orientation. For a given chem- 
ical potential of the passivating species the preference of a cer- 
tain rod type roughly follows the trend that the surface energy 
of the wire crystal is lowered by decreasing the density of sur- 
face dangling bonds. However, for extreme preparation con- 
ditions, e.g. in the presence of reactive atoms this route can 
be violated and passivated structures with complete adsorbate 
overlayers are preferred. 

In conclusion, we have shown that the rod energetics near 
the thermodynamic equilibrium is ruled by their surfaces. The 
formation energies of those rods approach the energies of the 
corresponding free surfaces for large diameters. The orien- 
tation of the surface facets and the passivation of their dan- 
gling bonds play an important role, in contrast to the bilayer 
stacking in growth direction, which only leads to minor en- 
ergy variations. Our studies are able to explain qualitatively 
a wide variety of results concerning bilayer stacking and rod 
facets found for III-V semiconductors. At least clear trends 
can be derived for thermodynamic equilibrium conditions. In 
addition, the growth kinetics and the catalyst, the growth con- 
ditions and after-growth treatment, e.g. annealing and the cor- 
responding atmosphere, may have an influence on the result- 
ing nanorod. 
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